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iy~j ' Kraichnan's model of passive scalar advection in which the driving velocity 

field has fast temporal decorrelation is studied as a case model for under- 
standing the appearance of anomalous scaling in turbulent systems. We 

o 

' demonstrate how the techniques of renormalized perturbation theory lead 

^ ' (after exact resummations) to equations for the statistical quantities that 

>' 

reveal also non perturbative effects. It is shown that ultraviolet divergences 
in the diagrammatic expansion translate into anomalous scaling with the 
inner length acting as the renormalization scale. In this paper we compute 
analytically the infinite set of anomalous exponents that stem from the ul- 
traviolet divergences. Notwithstanding, non-perturbative effects furnish a 
possibility of anomalous scaling based on the outer renormalization scale. 
The mechanism for this intricate behavior is examined and explained in de- 
tail. We show that in the language of L'vov, Procaccia and Fairhali [Phys. 
Rev. E 50, 4684 (1994)] the problem is "critical" i.e. the anomalous ex- 
ponent of the scalar primary field A = A c . This is precisely the condition 
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that allows for anomalous scaling in the structure functions as well, and we 
prove that this anomaly must be based on the outer renormalization scale. 
Finally, we derive the scaling laws that were proposed by Kraichnan for 
this problem, and show that his scaling exponents are consistent with our 
theory. 

PACS numbers 47.27.Gs, 47.27.Jv, 05.40.+j 
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1. INTRODUCTION 



The model of passive scalar advection with rapidly decorrelating velocity field 
which was introduced some time ago by R.H. Kraichnan M was suggested recently 

as a case model for understanding multiscaling in the statistical description of 
nonlinear field theories. The model is for a scalar field 0(r, t) where r is a point 
in R d . This field satisfies the equation of motion 

[d t - «V 2 + u(r,t) • V]6(r,t) = f(r,t). (1.1) 

In this equation /(r, t) is the forcing and u(r, t) is the velocity field which is taken 
to be a stochastic field, rapidly varying in time. The forcing is taken to be delta 
correlated in time and space-homogeneous, 

(f(v,t)f(r',t')) = %(v-v')5(t-f). (1.2) 

We will study this model in the limit of large Peclet (Pe) number, which is defined 
as the dimensionless ratio UlL/k, where Ul is the scale of the velocity fluctuations 
on the outer scale L of the system. The most important property of the driving 
velocity field from the point of view of the scaling properties of the passive scalar 
is the tensor "eddy diffusivity" [[|] 

hij(B.)= / dr([ui(r + K,t + T)-Ui(r,t + T)] 

J 

x [ % (r + R,t)- M ,M)]>) (1.3) 

where the symbol (...) in Eq. (|1.3j ) stands for an ensemble average with respect 
to the statistics of u which is given a priori. The scaling properties of the scalar 
depend sensitively on the scaling exponent (h that characterizes the R dependence 



of hij(H) 



h ij (R) = h(R) (a + i)^-a A;/ ' y 



R 2 



:i-4) 



h(R) = H[-) , (1.5) 



where C is some characteristic scale of the driving velocity field. The structure of 
the tensor is dictated by the incompressibility of the velocity field, and H is 
a free dimensional parameter. By "scaling properties" we mean how various cor- 
relation and response functions of 0(r, t) and its gradients depend on separation 
distances. For example the structure functions of 0(r,t) are 

^ n (R) = ([0(r + R,t)-e(r,t)] 2 "') . (1.6) 

In writing this equation we assumed that the statistics of the velocity field leads to 
a stationary and space homogeneous ensemble of the scalar 0. If the statistics is 
also isotropic, then S 2n becomes a function of R only, independent of the direction 
of R. The scaling exponents of the structure functions S2 n {R) characterize their 
R dependence in the limit of large Pe, 

S 2n (R)~R <:2 " } (1.7) 

when R is in the "inertial" interval of scales that will be discussed later in this 
paper. 

In Ref. Kraichnan showed that when the driving velocity field u(r, t) is 
delta correlated in time one can derive the exact "balance" equations for the 
structure functions S 2n (R): 

- B(R)S 2n (R) = J 2n (R) . (1.8) 

In this equation B(R) is the linear operator that will be used below repeatedly: 

B(R) = -R 1 - d ^R d - 1 h(R)^, (1.9) 

with d being the space dimension. On the RHS of the balance equations we have 

J 2n (i?)=4nK([0(r + R)-0(r)] 2 "- 2 |V0(r)| 2 ) . (1.10) 

Kraichnan conjectured that the scaling dependence of J 2n {R) on R when R is in 
the inertial range is given by the law 
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J2n(R) = J 2 {R)S 2n {R)/S 2 {R) ~ R^-b . ( L11) 

This scaling law led Kraichnan to far reaching conclusions regarding the scaling 
exponents ( 2n . Once inserted in the balance equations this conjecture resulted in 
the prediction that the scaling exponents satisfy the equation 

C2n(C2n-C2 + d) = ™*C2 . (1.12) 

If so, this model may be the first nonlinear non-equilibrium case where "mul- 
tiscaling" can be explicitly demonstrated. Further, Kraichnan et al. proposed 



in || that if this model is multiscaling, then ( |1.12| ) is the unique solution for 
( 2n . We will see that among other things our considerations lead to precisely the 
scaling law J 2n (R) ~ R& n ~& as conjectured by Kraichnan. We will see how the 
possibility of multiscaling is indeed realised due to the analytic structure of the 
theory. 

In a previous work on this model 0] (referred to hereafter as paper I) renor- 
malized perturbation theory was employed to study the scaling behavior. It is 
appropriate to present first a short summary of the results of this paper. 

The main point of paper I was that various statistical quantities exhibit 
anomalous exponents that stem from ultraviolet divergences in their diagram- 
matic expansion. For example there exists a quantity called the nonlinear Green's 
function (that is defined in sec. 2 and is a 4-point quantity) and which is denoted 
as {?2(0|i"i, r 2 , r3, r-4). It was shown that for r± ~ r 2 — r, ~ ~ R, and R^> r 
(cf. I, Eqs.(4.11), (4.12) and (4.18)) 

d d 

a a </2(0|ri, r 2 ,r 3 ,r 4 ) ~ r~ A , (1.13) 
ori or 2 

where A is an anomalous exponent that characterizes the leading divergence when 
r —>■ 0. Moreover, the value of this exponent is important in determining much 
of the scaling behavior in this model. It appears prominently in J 2n and also in 
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the correlation function of passive dissipation fluctuations. The dissipation field 
is defined here as 



e(x) = k|V6(x)| 2 , (1.14) 

and the correlation function K(R) is 

K(R) = (e(r + R,t)e(r,t)) - (e(x)) 2 . (1.15) 

It was shown in paper I that the ultraviolet divergences resulted in a dependence 
on the inner cutoff of the theory, denoted as rj [and defined below in Eq. ( |2.25|) 1 
which is written as: 

K (R)/k 2 ~ tT 2A (1.16) 

It was also explained in paper I that the theory indicates that if A reaches the 
critical value A c = special considerations must be made. For A < ^ it was 
shown that the perturbation theory for S4 and higher order structure functions 
converges order by order both in the ultra-violet and in the infra-red limits. 
Accordingly, there is no external renormalization length scale in the theory and 
there is no (perturbative) mechanism to ruin simple scaling, i.e. C,2n = n( n . In 
this case the correlation of dissipation fluctuations can be shown to decay in the 
inertial range of scales as 

( 20, -2 A 



K(R) ~ <e(*)> 2 ('!') (1.17) 



Indeed, as long as A < (h the correlation decays in the inertial range as it must for 
a mixing field. On the contrary, if A = (h Eq. (|1.17 ) cannot continue to hold since 



it predicts that the correlations do not decay. This is precisely where the need 
for anomalous scaling of the structure functions comes in. We will see below that 
the value of A is precisely (h, and we will get instead of Eq. ( |1 . 1 7| ) the following 
prediction: 



K(R) ~ (e(x)) 2 (|) 2 4 . (1.18) 

Here L is the outer renormalization scale for this model, that in general is not 
identical with the outer scale C of Eq.( |1.5D . This result is in agreement with 
Kraichnan's conjectures. The detailed explanation of how this phenomenon takes 
place is one of the major aims of this paper. It is an important mechanism 
that indicates how at least in this example ultraviolet divergences coupled with 
nonperturbative effects may conspire together to give at the end anomalous cor- 
rections which are carried by the outer renormalization scale. Whether or not 
such a mechanism operates in other hydrodynamic systems will be discussed in 
separate publications. At any rate, the example treated here serves as a powerful 
demonstration of the fact that renormalized perturbation theory as applied to 
field theories of the hydrodynamic type allows one, after proper resummations, 
to capture subtle nonperturbative effects. 

It needs to be stressed that the anomalous exponent A discussed above is 
just the leading divergence associated with scalar anomalous fields. We will ex- 
plain below that there exists a full spectrum of anomalous exponents which are 
associated with the inner length, and they have to do with anomalous fields of 
different irreducible representation of the rotation group 0. We will compute be- 
low analytically the whole spectrum of these exponents. For the model at hand 
this is an easy task, but it serves to demonstrate the rich scaling properties of 
hydrodynamic systems. This rich scaling structure has not been considered by 
the fluid mechanics community until now. 

The structure of this paper is as follows: section 2 is devoted to the derivation 
of the differential equations satisfied by the n-point time correlations functions 
and the n-point simultaneous correlation functions, as well as the equations for 
the 2-point and 4-point Green's functions. The derivation is based on the di- 
agrammatic expansion of paper I, and the main objective of this section is to 
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demonstrate that this technique yields the exact equations, and that the result- 
ing equations contain aspects of the problem which are not perturbative. The 
results of the calculations of this section are identical to alternative derivations 
which are based on standard stochastic methods. Readers who are not interested 
in the method of derivation can start to read this paper from section 3 which 
begins with a catalog of the equations that are analyzed in the rest of the paper. 
Section 3 begins with the exact solution of the Green's function and the 2-point 
correlation function. These solutions are important since they introduce the ho- 
mogeneous solutions of the operator B{R) which then appears importantly in all 
the solutions of the higher order quantities. Section 2 B deals with the exact 
solution of the 2-point correlation function. This solution depends on the nature 
of the forcing. We are able to study in detail how the effects of nonisotropic 
forcing on the large scales decay in the inertial range as the scale of observation 
decreases. The exponents that govern the "law of isotropization" are important 
in determining the scaling behavior of correlation functions of anomalous fields 
whose presentation under the symmetry groups is different from scalars. Finally, 
in section 2 C we discuss the 4-point nonlinear Green's function which allows us 
the evaluation of the anomalous exponent A. Section 4 is devoted to the calcu- 
lation of K(R), J±{R), and other correlations that expose non-scalar anomalous 
fields. The section is based on analyzing the equation for JF 4 , and the strategy 
is to extract the leading divergence that is characterized by the anomalous ex- 
ponent A. Section 5 presents a brief calculation of J 2n and section 6 collects 
all the result together in order to compute the scaling exponents of S , 2 „(i?). We 
show that there are only two possibilities, one is simple scaling, and the other 
is anomalous scaling in agreement with Kraichnan's conjectures. Simple scaling 
is possible only if the dissipation field is not mixing, which seems a nonphysical 
condition. Section 7 offers a summary of the paper. 
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2. DIAGRAMMATIC DERIVATION OF THE EQUATIONS FOR 
CORRELATION AND RESPONSE FUNCTIONS 



In this chapter we demonstrate that for this model the exact resummation 
of the diagrammatic expansion for the various n-point correlation and response 
functions of the theory results in exact differential equations that contain both 
perturbative and nonperturbative anomalies. The analysis is based on paper I, 
in which the first step is the Belinicher-L'vov transformation [§-0]. This is done 
by allowing the center of the coordinate system to move along the Lagrangian 
trajectory of a particular fluid point. The reference point is at position ro at time 
to- The trajectory of this point with respect to r is 

p(r , t \t) = f dt' u(r + p(r , t \t'),t'). (2.1) 

J to 

Let us denote the transformed variables as 

T(r , t \r, t) = 0(r + p(r , t \t), t), (2.2) 
v(r , t \r, t) = u(r + p(r , t \t), t), (2.3) 
0(r o , t \r, t) = /(r + p(r , t |*),*)- (2.4) 

In terms of these new variables the equation of motion ( |1 . 1|) reads 

T(r , t \r, t) = 0(r o , t \r, t), (2.5) 

where 

w(r , t |r, t) = v(r , t |r, t) - v(r , t \r , t). (2.6) 

We turn now to the discussion of the statistical quantities that are defined in 
terms of these variables. 

A. The 2-point quantities 

The 2-point functions discussed here are the 2-point correlation function 
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d t - kV~ -w(r ,t |r,t)-V 



J r (r \x 1 ,x 2 ) = {T(x q \x 1 )T(xq\x 2 )) (2.7) 



and the Green's function that is defined as 



G{r Q \x 1 ,x 2 ) 



/ 5T{x (y \x l ) 



(2.8) 



\ 5£(x \x 2 

where the d + 1 dimensional vector x = (r,t), and £ is the Belinicher-L'vov 



transformation of a forcing that is added to the RHS of Eq. (|2.5| ). The dependence 
on to disappeared from the LHS of Eqs.( p.8|) because one can prove 0] that all 
the average quantities are time translationally invariant. 

In this subsection we present the derivation of the equations for ^(rol^i, x 2 ) 
and Q(r \xi, x 2 ). In paper I we derived the Dyson- Wyld equations for these two 
point functions. For t > they read 

(d t - /cV 2 )0(r o |ri, r 2 , t) = 5{y x - r 2 )5(t) 

+ J dr , V(T Q \r 1 ,i')G(T \r f ,r 2 ,t) (2.9) 

/roc 
dx'dx" J dt'g(ro\r u r',t + t') 

x [%(r',r") + ^r',r")]g(r \r 2 ,T",t'). (2.10) 

For negative times the Green's function is zero, and the correlation function is 
symmetric to inverting time t — > — t and the coordinates r — ► — r together. In the 
Dyson equation ( |2.9|) the mass operator S(ro|r,r') can be written explicitly 

E(ro|r,r , ) = ^(ro|r,r / ,t = 0)^(r |r,r / )^ 7 (2.11) 



whereas in the Wyld equation (|2.10| ) the mass operator $ takes the form 



$(r |r,r') = ^(r |r,r')^J 7 ^(r |r,r / ,t = 0). (2.12) 



In flrilD and (|2TT2D 



/oo 
dt(wi(r , t \r, t)wj(r , t \r', 0)). (2.13) 
-oo 
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In the case considered in this paper in which the velocity field is delta correlated 
in time the expression ( |2.13| ) can be related to the eddy diffusivity ( |1.3| ), expressed 
in terms of Eulerian correlations as: 



Hij(r \T, r') = hij(r - r ) + h^r' - r ) - /^(r - r'). 



(2.14) 



That the mass operators (|2.11|) and (|2.12| ) have an explicit form corresponding to 
the 1-loop diagram rather than an infinite series representation stems from the 
fact that the velocity field decorrelates on an infinitely short time scale. This 
simple form of the mass operators will be lost if we relax this fast decay of the 
velocity time correlation functions. 

Equations ( |2.9|) and ( |2.10| ) can be turned into differential equations for the 
2-point functions T and Q. Substituting ( [2.11| ) 



in 



we find that 



(d t - KV 2 )£(r |r,r',t) = 5(r - r')5(t) 



(2.15) 



+ J dr"-^g(r \v, r", 0)^(r |r, r ")^G(r Q |r", r', t). 







Remember that Q(ro\r, r', t) is zero for negative times. Integrating ( [2.15 ) over 
time from t = —5 to t = 5 and taking the limit 5 — > one finds that 



£(r |r,r',t = + ) = 5(r - r') . 



(2.16) 



We will choose symmetrical regularization at t — and by convention write 



g(T \T,T>,t = 0) = -6{T-I>). 



(2.17) 



Using this evaluation the integration may be performed leading to 



d t +A(r- r )l G(r \r, r', t) = 5(r - r')5(t). 



(2.18) 



Here we introduced the generalized diffusion operator 



^(r) = -kV 2 + B(t) 



(2.19) 
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where the operator B is a key operator that appears repeatedly below. We will 
distinguish between an operator B(r a , r^) which acts on functions of two variables 
r a and and an operator B(R) acting on functions of one variable R: 

d 2 

B(r a , Tp) = B a> p = h i:j (r a - 173) — , (2.20) 

Clearly B(r a , rg) turns into £>(R) on the class of functions depending on the 
difference R = r a — rp only. In spherical coordinates the S-operator can be 
represented as the sum of two contributions: 

B(R) = B{R) + (C " ( ;' 1) 1) ^ • (2-22) 



Here B is the operator ( |1.9| ) and L is the angular momentum operator —ir x V 
which depends only on the direction of r. 

For future purposes we also need the equation of the time-integrated Green's 
function Q(r — r , r' — r ) which is defined as 

0(r-r o ,r'-r o ) = J dtG(r \v,r',t). (2.23) 

This function satisfies the equation 

T>! (R)^(R, R') = S(R - R') , (2.24) 

which follows from Eq. (|2.18|) . 

The equation of motion (|2.24|) allows us to introduce the inner scale of this 
model, denoted by rj. By definition rj is the scale for which the advective term B 
is of the order of the dissipative term kW 2 . Equating the two terms we get 

/ K \ 

V-[jj) ■ (2.25) 

We proceed now to determine an equation of motion for the two-point corre- 
lator .F(ro|ri, T2, t). It is clear from Eq. ( |2.18p that one may define an inverse 



operator for the Green's function C/(ro|r, r', t) according to 
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g^(r-r ,t) = d t -A(r-r ) . (2.26) 

Note from Eq. ( |2.18|) that the equation of motion for Q(r \r,r',t) only depends 
on the first coordinate r, so that (Jf 1 = Q^ij — io,£). Operating with Q^ 1 we 
may rewrite the Wyld equation (|2.10|) as 



d t + T>i(r 1 -r )] [-<9t + X>i(r 2 - r )] J r {r \ri,r 2 ,t) 
= 5(t) [$(r |r 1; r 2 ) + $ (r |ri, r 2 )] , (2.27) 

where we have used the fact that jF(ro|ri, r 2 , t) is only a function of the time 
difference t — t\ — t 2 . We have also the "boundary" condition 

^ 7 (r |r 1 ,r 2 ,*) -»• for t -> ±oo . (2.28) 

Next we want to derive the differential equation satisfied by the simultaneous 2 
point correlator JF(R): 

^(R) = ^( ri ,r 2 ,t = 0) (2.29) 

where R = ri — r 2 . The derivation is described in detail in Appendix A with the 
final result 

X> 2 (R).F(R) = $ (R), (2.30) 

where generally the operator ^(r,,,^) operates on two coordinates: 

£> 2 (r Q , r/s) = -K[V 2 a + Vj] + B(r a , v p ) , (2.31) 

but in the case where the operand is a function only of the difference R = r Q — r@ 
reduces to 

£> 2 (R) = -2kV 2 + B(R) . (2.32) 
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B. The derivation of the differential equations for higher order correlations and 

response functions 



In this subsection we present the equations for the 4-point and higher order 
correlation functions and for the nonlinear Green's function ( |2.8| ). The solutions 
are deferred to the next section. 



1. The 4-point Green's function 

In paper I we presented the diagrammatic series for the nonlinear Green's 
function Q 2 (r \xi, x 2 , x 3 , x 4 ). This quantity is defines as 



£ 2 (i"oki, x 2 , x 3 , xa) 



/ 8T(x \xi) 5T(x \x 2 ) \ 



(2.33) 



\ S^(x \x 3 ) $£(x \x 4 ) I 
The diagrammatic expansion of this quantity is an infinite series of ladder 
diagrams that can be resummed exactly. In Fig. 1 we recall the notation for the 
diagrammatic elements, and in Fig. 2 is reproduced the diagrammatic resummed 
equation for this function. In analytic form this equation reads 

£2(0|£i, x 2 , £3, x 4 ) = Q 2 (0\xi,x 2 ,x 3 ,x 4 ) 
+ f dr'dv" r dt' g 2 {0\ Xl ,x 2 , r', t', r", t') (2.34) 

J J £ m 

xHdv', 0^7^(01^ t',r",t',x 3 , x 4 ) 

where t m = min{t 1 ,t 2 } and 

Q 2 (0\x u x 2 , x 3} x 4 ) = g(0\x 1} x 3 )g(0\x 2 ,x A ). (2.35) 

We want now to derive a differential equation for this quantity. We can proceed 
as in the case of the correlation function by applying the product of the inverse 
operators 1 of Eq. ( |2.26| ) to obtain a differential equation for Q 2 {$S\x\, x 2 ,x 3 , 24) 
as a function of three time differences. We may however make use of the structure 
of Eq. (|2.34j) to directly derive an equation in only one time difference. The 
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availability of such a reduction follows from the delta-correlated velocity field in 
much the same way as the availability of a differential operator for .F(0|ri, r 2 , t = 
0) as we saw before. We can choose at will to consider Eq. ( |2.34| ) for the times 
t\ = ti = t and £ 3 = £4 = to derive a differential equation in one time difference 
for the quantity 

£ 2 (0|r 1; r 2 , r 3 , r 4 , £) = £ 2 (0|r 1; £, r 2 , £, r 3 , 0, r 4 , 0) (2.36) 

The derivation resembles closely the derivation of the equation for the two-point 
correlator described in Appendix A. Applying the operator ( |2.27| ) to Eq. ( |2.34p 
with the above choice of times we get 

[d t + P 2 (n, r 2 ) + H{t u r 2 )]^ 2 (0|n, r 2 , r 3 , r 4 , £) (2.37) 
= 5(ri - r 3 )5(r 2 - r 4 )5(£) , 



where £> 2 (i"i, r 2 ) is the operator defined in Eq. ( |2.32| ) and 7i is an operator defined 



in Appendix A that arises from the derivation of the equation for the two-point 
correlator. 

For the time integrated 4-point Green's function 

£ 2 (r 1; r 2 , r 3 , r 4 ) = J dtg 2 (r 1 , r 2 , r 3 , r 4 , £) (2.38) 

one has the equation 

X> 2 (n, r 2 ) + W(ri, r 2 )] &(n, r 2 , r 3 , r 4 ) = 6(n - r 3 )6(r 2 - r 4 ) . (2.39) 

The non-linear Greens' function is the kernel of the response to some forcing 
at points r 3 and r 4 , and we are interested in this response. In order to study this, 
we introduce a new function 

(r la r 2 ) = J dr 3 dr 4 g(r u r 2 , r 3 , r 4 )A(r 3 , r 4 ) , (2.40) 

where A(r 3 , r 4 ) is an arbitrary function. From Eq. ([2.39|) one may determine that 
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P 2 ( ri ,r 2 )+ft *(n,r 2 ) = A(r 1; r 2 ). 



(2.41) 



If one chooses now to restrict A(ri,r 2 ) to the space of functions depending only 
upon the difference R = ri — r 2 , the equation simplifies (as in the derivation in 
Appendix A) to the form 



P 2 (R)*(R) = A(R). 



(2.42) 



We will return to an analysis of this equation in Sec.3C. 



2. The simultaneous higher order correlations 

The equations of motion for the time dependent 2n-th order correlation func- 
tions are derived in Appendix B. Here we will derive the equations for the simul- 
taneous correlations. The simultaneous 2n-point correlator is 

F 2n (n, r 2 , r 2n ) = (T(0|ri, t)T(0\r 2 , t), T(0|r 2n , t)) . (2.43) 

This same time quantity is identical in the Eulerian frame of reference and in the 
transformed reference frame that we use. We can thus forget the r designation, 
and remember that in homogeneous systems the quantity is a function only of 
differences of its space arguments. Its time derivative is on the one hand zero and 
on the other hand 

d 

^2n(ri,r 2 ,...,r 2 „) (2.44) 

2n I d \ 

= WT(0| ri ,t) ... ^T(0|r Q ,t) ... T(0\r 2n ,t)j . 

Using the equation of motion ( |2.5|) we find that the RHS of (|2.44|) has three 
types of terms, one with w (advection), one with k (dissipation) and the last 
proportional to the force (forcing). These terms are denoted as A^v, ^4dis and 
A ior with 
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+ ^-dis + -^for — 



(2.45) 



where 



2n 



^adv — E ^ 7 a ' ^w,2nT{ r a, r l--- r 2n) 



(2.46) 



a=l 



and 



(w a T(0| rij t) . . . r(0|r a _ ls t)T(0|r Q , t) . . . T(0|r 2n , t)) 

2n 

A dis = « E V a^2r.(ri, r 2 , ... r 2n ) . 



(2.47) 



(2.48) 



Finally 



2/i 



^for = E <^(0|n, t)...T(0|r a _i, t)0(O| rQ5 t)...T(0|r 2n , t)) . 



(2.49) 



The diagrammatic representation of F Wj2n T( r a, r i--- r 2n.) wa s discussed in paper I 
with a final result [cf. paper I Eq. (5.7)] which in r, t representation is 



F w OnT\^a, ri---I"2r. 



(2.50) 



2/i 



- H (0|r a , Tp) ■ V /3 J r 2n(ri, r 2 , r 2n ) 



Introducing this result into ( [2.47[ ) and using (|2.14 ) we find 



2n 



(J 2 



X 



-J~2n( r li r 2, I"2n)- 



draidrpj 

This form can be rewritten equivalently as 



(2.51) 



A 



adv 



2n 

E 4/3 

a>P=l 



[2.52) 



+ 



2// 



Eh 



o 



La=l Ul m 



2n o 
V — 



^2n( r lj r 2, *2n) 



Finally we use the fact that Tm is a function of differences of its spatial arguments 
to see that the second operator in the first line of the RHS gives zero and the line 
can be omitted. Thus 
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2/i 



^4adv — — X] &ctpJ~2n{Yl) *"2, r 2n )- 

a>P=l 

Consider next the forcing term. Using the fact that for a Gaussian force 



(T(0|n, *)...T(0|r tt _i, t)0(O|r Q , t)...T(0\r 2n , t)) 
^T(0|n, t)...T(0|r a -i,t)...T(0|r 2n , t) \ 



(2.53) 



x $ (r Q - r^)rfr^ 



(2.54) 



where we used (1.2) and the fact that we deal with simultaneous correlations. 
The functional derivative in the integrand is a zero time response, which as usual 
is computed in the non-interacting limit: 



/ 5T{0\t u t)...T(0|r Q _ 1 , t)...T(0|r 2ra , t) \ 

\ WIC *0 / 

^ G°(r Q , — Tp, t = 0)JF 2n~2( r l, r 2, ••• r 2n) 



(2.55) 



where in J-2n-2( r i, i"2, ■■■*in) the two arguments r a and vp are missing. Substitut- 
ing ( |2.55p in Q2.54|) and the result in ( |2.49| ), we find 



A {oi = $ o(r Q - i7j).F 2n _ 2 (ri, r 2 , ...r 2n ). 

a>/3 



(2.56) 



Substituting all these results in ( |2.45|) yields 



2 1 1 



•7"2n( r l> r 2, •••) I"2n 



a a>/3 

= ^ $o(r a - r^) ^2n- 2 (ri, r 2 , ...r 2n ) • 

a>,3 



(2.57) 



This equation, which is our final equation for the simultaneous correlation func- 
tion, is again identical with Kraichnan's. 

3. Equation for the irreducible correlation T% 
The irreducible (cumulant) part of the correlation functions are defined as 



*77C rp pG 

J ~2n ~ ^In J ~2n- 



(2.58) 
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The equation for the simultaneous T% follows from specializing Eq. ( |2.57 ) to the 
case 2n = 4, using Eq. ( |2.30| ) for JF 2 , 



E«^+E Bap 

a a>/3 



J2 ^4( r b r 2,r 3 ,r 4 ) 

a>/3=l 



= -« E - r a3 )J r 2 (ri,r 2 )J'2(r3,r4) 

(a.) 

pcrm(1234) 

= - '(B14 + #23 + Big + 524)^12^34 
+ (Bl4 + #23 + #12 + #34)^13^24 
+ (Bl3 + #24 + B 12 + 44)^14^23" • 



(2.59) 



3. ANALYSIS OF THE GREEN'S FUNCTIONS AND THE 2-POINT 

CORRELATION 

In this section we begin to discuss the exact solution for the various quantities 
in this model when all the separation distances are in the inertial interval. For 
the sake of clarity we catalogue here all the equations that we are going to solve 
below, neglecting the diffusive terms which appear in the full equations. The 
equations are the following, for 

1) the time integrated Green's function cf. ( |2.23| ): 

5(R)£(R, R') = 5(R - R) , (3.1) 

which follows from ( [2.24 ), 

2) the 2 point simultaneous correlation functions, 

B(R).F(R) = $(R) , (3.2) 

which follows from (|2.30 ), 

3) the function \I/ introduced in Sec.2B, 

B(R)*(R) = A{R), (3.3) 
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which is the k — > limit of Eq. ( 2.42|) , 

4) the cumulant of the 4-point simultaneous correlation function 

B12 + Bja + B u + B 23 + B 24 + i3 34 ]^ c (ri, r 2 , r 3 , r 4 ) 

B u + B 23 + Big + B 24 ]F{ ri - r 2 )^(r 3 - r 4 ) 
B u + B23 + B 12 + B u ]F( ri - r 3 )^(r 2 - r 4 ) 
Bia + S 24 + B12 + i3 3 4l^(ri - r 4 )^(r 2 - r 3 ) , 



(3.4) 



which is the k — > limit of eq. ( [2.5S ) . 

It is noteworthy that the same operator B, (|2.20|) , ( p.22j ) appears in all these 
equations. Indeed, the most interesting scaling properties will be seen to arise 
from the eigenfunctions of B with zero eigenvalue which are the solutions of the 
homogeneous equation (j3.5|). 



A. Solutions of the basic homogeneous equation 

For two distinct coordinates R ^ R' the Green's function Q (R, R') satisfies 
the homogeneous part of the equation ( p.l[) 

(3.5) 



£(R)£(R,R') = , 



which will be referred to hereafter as the " basic homogeneous equation" . In light 
of the representation of Eq. (|2.22j ) of the B operator via the angular momentum 
operator we represent the solution of (|3.5|) in 3-dimensional space as an expansion 
over spherical harmonics 



oo I oo I 

<7(r,r') = £ E E E 

"ll',mm' 9lm 

1=0 m =-ll'=0 m'=-V 

where bw mm i are coefficients and 



(3.6) 



9h 



(R) = R^Yi 



lm \ 



(3.7) 
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The spherical harmonics Yi m (9,4>) are the eigenfunctions of the angular momen- 
tum operator: 

L 2 Y lm (9,<f>) = l(l + l)Y lm (e,<f>) (3.8) 

and 9 and 4> are the polar and azimuthal angles of R. Substituting the expansion 
])into eq. (|3.5|) we find the relationship 



Wi + 3 - Ca) - Kl + 1)(1 + a/2) = . (3.9) 

Of the two solutions of this quadratic equation we must select the nonnegative 
branch since the negative branch is unphysical: 



C 2 -3 + V(C2-3) 2 + 2/(/ + l)(2 + a) 



(3.10) 



In 2 dimensions one has a similar representation with the difference that instead 
of the spherical harmonics we have the eigenfunctions exp(il(p) and instead of the 
eigenvalues 1(1 + 1) we have the eigenvalues I 2 . The final result for the exponents 
Pi in 2- dimensions is 



a + va + 4^(1 + go 



(3.11) 



Multipole expansions of the type ( ^.6| ) and the scaling exponents $ play an 
important role in all our development below. We will see that the basic homo- 
geneous equation reappears in various guises, and in each of them we will need 
to fix the coefficients bu' tmm t such as to respect the symmetry of the particu- 
lar object involved, the boundary conditions, etc. For example in this case we 
know that the Green's function C?(R, R') is symmetric in R and R'. Accordingly 
bii',mm' = bi'i,m'm- If the solution depends on R, R' and the angle between R and 
R', then bw mm / ^ only for 1 = 1', m + m! = 0, etc. 



B. Exact solution of the 2-point correlation function 
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1. Isotropic forcing 



The solution of the 2-point correlation function depends on the nature of 
the external forcing f(r,t). It is customary to take f(v,t) to be Gaussian and 
statistically homogeneous in space and time. The properties of the correlation 
function $o( r — r ') were not determined up to now. Since we are interested in the 
universal scaling properties of the scalar field we want to choose the forcing such 
that it has only large scale components. Otherwise the scaling exponent of the 
2-point correlation may be coloured by the functional dependence of the forcing 
on r. On the other hand the forcing may be isotropic or non isotropic. In this 
subsection we will deal with the isotropic case, and the nonistropic case will be 
treated in the next subsection. 

The properties of f(r,t) are best stated in k-space: it is concentrated in the 
small k region, i.e. k < 1/L, and it decays quickly to zero for fc < 1/L. In r 
space this means that (/(r, t)f(r + R, t)) is constant for R <C L: 



Using this form of the forcing correlation function we realize that the inhomoge- 
neous term on the RHS of Eq. (|2.29j) is translationally and rotationally invariant. 
Accordingly we can seek solutions that have the same symmetry, solving the 
equation 



(/(r,t)/(r + R,t)) = <D . (R « L) 



(3.12) 



B(R)T{R) = $, 



(3.13) 



The solution of the inhomogeneous equation is 



J~inh(R) — CmhR'" 2 , 



(3.14) 



with C\ 



inh 



$ /2/<2 and 



C2 = 2 - Ch • 



(3.15) 
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Of course, the inhomogeneous solution has to be supplemented with the solutions 
of the homogeneous equation in order to match the boundary conditions (which in 
this case are ^(R) = ^"(O) 7^ at r = 0). There are two homogeneous solutions: 
one is a constant which must be taken as ^(O), and the other is F{R) oc i? - ^ 1 
which is ruled out by the boundary condition J^(0) < 00. One should note that 
the constant homogeneous solution belongs to the family of exponents (3i with 
/ = as one would expect. The solution that is ruled out is indeed the / = 
member of the forbidden branch of solutions of the quadratic equation (|3.9| ). 



2. The law of isotropization 



When the forcing is nonisotropic on the large scales we need to apply the full 
operator i3(R). In the regime R -C L the forcing is again a constant which is 
independent of the angles. However, for R ~ L we expect nonisotropic forcing 
which leads to nonisotropic tails for every value of R. We can find the law of 
isotropization of ^(R) for R small by solving the equation that involves now also 
the nonisotropic part of the operator B(R) 



B(R).F(R) = $ , R < L. 



(3.16) 



The inhomogeneous solutions are the same as before, but in order to match with 
the anisotropic behaviour at R ~ L we need to invoke the whole set of homo- 
geneous solutions gi m (R), Eq. (|3.7p . The solution in 3 dimensions is written 

as 

OO l / J^\Pl-C,2 



J-(R) = ^(0) + R^ C m + J2 E a im T 

l=2m=-l KLJ 



(3.17) 



To understand this result we note that all values of (3i are larger than ( 2 - in 
particular (for both d = 3 and d = 2) is 
1 



02 



(2-d+J(( 2 + d)z + 24( h 



(3.18) 
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For positive (h, f3 2 is larger than ( 2 . For larger /, f3i is even larger, and dependent 
on dimension. Thus all the anisotropic terms decay when Note that the 

coefficients ai m are nonuniversal and should be found by matching at R ~ L. The 
law of decay is however universal. 



3. Anisotropic structure functions 



The solution Q3.17D suggests the introduction of anisotropic structure functions 
via the definition 

S 2 ,im(R) = J d9d(j)Yi m (9, <fr)S 2 (R) . (3.19) 

Here I should be even due to the symmetry with respect to the inversion of R. 
These anisotropic structure functions display "clean" scaling behavior with the 
exponent fif 

S 2M (R)~S 2 (L)(^-Y' . (3.20) 



L 

We will see that the same scaling exponents feature prominently below 



C. The nonlinear Green's function and the anomalous exponent A 

In this subsection we discuss the nonlinear Green's function and the anomalous 
exponent A which is associated with it, cf. (|1.13|). To this aim we return to 



the equation for Eq.( |3.3| ). We realize that the solutions of this equation are 
identical to those discussed already in the context of the 2-point correlator. For 
a constant function A we can write the solution 

\&(R, R') = \l/(0) + C|R - R'| ?2 . (3.21) 

Taking the derivative with respect to R and R' as required by Eq. (|1.13|) we find 

(Vi ■ V 2 )tf(R,R/) ~ |R-R'| _Ch (3.22) 
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Comparing now with the definition of the anomalous exponent A in Eq. (|1.13|) we 
see that in the limit |R — R'| — > rj this quantity diverges like l/rj^ h and 

A = (h- (3.23) 

As explained before this value of A is the critical value A c , and the implications of 
this finding are explored below. We turn now to the appearance of the anomalous 
exponent in the 4-point correlator and related quantities. 

4. THE FOUR POINT CORRELATOR AND RELATED QUANTITIES: K(R), 

AND J 4 (R) AND L(R) 

In this section we present the analysis of the four point correlator 
J r (r 1 , r 2 , r 3 , r 4 ) and of the quantities that are related it, i.e the correlation of 
dissipation fluctuation K(R), Eq. (|1.15| ) and the RHS of the balance equation 



J 4 (-R), Eq. (|1.10| ). These last quantities are not exactly 4-point correlations, but 
they are obtained as a limit of J-a(ti, r 2 , r 3 , r 4 ). The simpler limit is K(R) which 
is a centred correlation function and is therefore related to ^(ri, r 2 , r 3 , r 4 ): 

K(R)=k 2 lim lim (Vi • V 2 )(V 3 • V 4 ) 

x >*(n, r 2 , r 3 , r 4 ) + ^ 2 (n, r 3 )^ 2 (r 2 , r 4 )l . (4.1) 



Since we discovered that the the second derivative is singular , cf. Eq.( |3.22| ), we 
need to carefully examine the above limit. 

The quantity J 4 (i?) has a Gaussian decomposition Jf(R) and a cumulant part 
Jl{R). The Gaussian decomposition is trivially computed as 

jf(R) = eS 2 (R) + ^[VS 2 (R)} 2 . (4.2) 

The cumulant part is related to T% via two terms, J%(R) = J 4j i(i?) + Ja,2(R), 

J 4il ( J R) = K ((|V 1 0(r 1 )| 2 6 2 (r 2 ))) , (4.3) 
J4, 2 (^) = -2/€((|V 1 0(r 1 )| 2 0(r 1 )0(r 2 ))) , (4.4) 

25 



where double brackets denote the cumulant part. These two contributions may 
be considered as the following limits of J rc (r 1 , r 2 , r 3 , r 4 ) 

Ja,i(R) = -k lim lim (Vi ■ V 2 )^(ri, r 2 , r 3 , r 4 ) , (4.5) 

r 12, r 34^0 n 3 -+i? 

J4,., OR) = 2k lim lim (Vi ■ V 2 ) J^(t u r 2 , r 3 , r 4 ) . (4.6) 

ri2,ri3-+0 ri4— >R 

Note that in Eq.( |4.5| ) we have two pairs of coalescing points, i.e 1,2 and 3,4, which 
are separated by a large distance R. In Eq.( |4.6| ) we have three coalescing points, 
i.e. 1,2,3 and this group is separated from point 4 by R. One should also note 
that J 4 is obtained from the full JF 4 and not from the cumulant part. 

In the next subsections we are going to make strong use of the divergence with 
respect to small distances. Our strategy will be to expose the leading exponent in 
the divergence with respect to small separations and to compute it exactly. Then 
we will find the exponent of the dependence on R by power counting, knowing 
the overall scaling exponent £ 4 of T%. In other words, our basic assumption is 
that the correlator T% is a homogeneous function of its arguments as long as all 
of them are in the inertial range: 

^(Ar 1; Ar 2 , Ar 3 , Ar 4 ) = X^T%(v x , r 2 , r 3 , r 4 ) , (4.7) 

where £4 is the unknown scaling exponent that characterizes the structure function 
£4 OR). We will not make any assumption about the numerical value of £4. 



A. Two coalescing pairs of points 

1. The effective equation for T\ 



Consider Eq. (|2.59| ) in the limit r 42 , r 34 <C R, but all separations in the inertial 



interval. This allows us to neglect the diffusion terms, and to conclude that 
oc r and i3 3 4 oc r 34 '' 2 are much larger than all the other S a p operators 
which are proportional to R~&. This suggests to rewrite the equation in the form 
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Bt2 + Bza + E ^4(ri, r 2 , r 3 , r 4 ) = RHS , (4.8) 

where the sum is on all the four B operators other than the ones explicitly dis- 
played. The calculation of the RHS of Eg. ( j2.59| ) is elementary since we know 
everything explicitly. The result of the calculation is 



RHS = R C2 



(4.9) 



where A , Ai and A 2 are some dimensionless constants. The last term in the 
bracket comes from the first group of terms in the RHS of Eq.( [2.59| ). Indeed, 



the product of two JF 2 s is proportional to (r 12 r 34 ) < > 2 . The naive evaluation of the 
sum of jB's is R^ h /r 12 r 34 . However there are two cancellations in the combination 
of £Ts which result in the replacement of r 12 r 34 by R 2 . The sum of i3's becomes 
proportional to R~&. The leading contribution in the last two groups of terms 
is regular in ri 2 and r 34 and therefore can be evaluated as R&. This is the first 
term in the bracket. The second term in the bracket is contributed by the sum of 
&i2 + ^34 in the second and third groups of terms. Higher order terms exist but 
they are proportional to the second power of the small distances. 

To discuss Eq. Q4.8|) further we note that in a space homogeneous situation T1 
is a function of six differences in 3 dimensions and 5 in 2 dimensions. In light of 
our strategy it is convenient to choose the variables as R, r 42 and r 34 . In doing 
so we are using seven rather than the needed number of variables, but we will see 
that the dependence on the extra angle variables disappears. We can now group 
the RHS of Eq.( [4.8| ) together with BT% into a new function, say i?(r 12 , r 34 , R) 

£(r 12 , r 34 , R) = RHS - £ BT A . (4.10) 

In the limits ri 2 , r 34 <C R we can expand i?(r 42 , r 34 , R) in orders of r 42 and r 34 

12, r 34 , R) + ... , (4.11) 

where 

27 



E^(R)= lim E(r 12 ,r u ,R). (4.12) 

ri2,r34^0 

The equation that we need to analyze takes on the form 



B 12 + #34 



J- 4 c (r 12 , r 34 , R) = EW(R) + £ (1) (r 12 , r 34 , R) + . . . • (4.13) 



and the explicit expression for ^-^(r^, r 34 , R) will be presented below. 

2. Solutions 

To leading order we have the effective equation: 

Si2 + 634] ^(ri 2 , r 34 , R) = E®\R) . (4.14) 
The leading inhomogeneous solution of Eq. fl4.14|) is found by inspection: 



^L(r 12 , r 34 , R) = Ci£<W(J2)(rg + rg) , (4.15) 

where C\ is a dimensionless constant. Using the overall scaling exponent (4 that 

was introduced in Eq. ( |4.7|) we can rewrite this solution as 

C2 C2 

*S*(ri2, r 34 , R) ~ S,(R) r -^±^ . (4.16) 

This inhomogeneous solution will be particularly important in the context of the 
calculation of J^{R). It does not contribute however to the evaluation of K(R), 
and for the the latter we need to find the next order inhomogeneous solution. The 
next order term on the RHS of Eq.( |4.13| ) stems from two sources. First is the A\ 
term in Eq. (fO|) , and the second is found by substituting ( |4.16| ) into J^RFI on 



the RHS of ( [4.121) . Both contributions have the same dependence on r\ 2 and r 3 4, 



and they can be written as: 

E^(r 12 ,r u ,R) = El (R) r -^±^ . (4.17) 

Substituting this in Eq. fl4.13p produces the next order inhomogeneous solution 
which reads 
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r 3 4, J2) = Ci^i?) (4-18) 



We can again use Eq. fl4.7f) to rewrite this solution in the form 

^L(ri2, r 34 , R) ~ 5 4 (i2) • (4-19) 

This solution will be shown to give the leading order contribution to K(R), and 
therefore these orders of the inhomogeneous solution will suffice for our analysis. 

In addition to the inhomogeneous solutions all the homogeneous solutions 
fl3l)| ) found below are available to us, since the homogeneous solutions of B\i and 
of S34 are also homogeneous solutions of Bu + $34. We can therefore write the 
homogeneous solution as 

00 1 00 1 

^4,hom(ri2,r 3 4,-R) = E E E E f ll',mm>(R)gim(r 12) gi>m> (r U ) . (4.20) 
1=0 m=-ll'=0 m'=-V 

One should note that at this point we can make a choice of the coordinate system 
such that the z axis is directed along the R axis. Accordingly the dependence on 
the sum of azimuthal angles 0i2 and ^34 should disappear due to the symmetry 
of the problem. This requirement is met if all coefficients fiv^mm 1 vanish when 
m + ml 7^ 0. In addition our correlation function is symmetric with respect to 
the exchange of any pair of points. Accordingly all odd values of I and V are 
excluded from the sums in ( OO) . Finally we can find the R dependence of these 



coefficients from the overall power counting. Since ^ m (r) oc r@ l we can write 

fw, mm '{R) ^S^R)/R^ +p i . (4.21) 

We will need below also the next order inhomogeneous solution which is ob- 
tained by substituting ( 4.20Q back into the definition of E ( f4.10|) . This solution 



will be denoted as 



C2 C2 

Fi% h (r 12 , r 34 , R) ~ *gL(ri2, r 34 , R)"- 2 ^ 1 ■ (4-22) 
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In addition to these homogeneous solutions one can consider also the homogeneous 
solutions of the sum of the two B operators on the LHS of ( |4.14| ). It can be 
shown that these homogeneous solutions do not add any information that is not 
contained in the solutions described above. 



B. Three coalescing points 

For the calculation of J4 we need also to consider the geometry of three coa- 
lescing points, see Eq.( |4.4|) . Accordingly we focus on the limit r 12 , r 13 , r 23 <C R ~ 
ri4 ~ r 24 — r 34 . Instead of Eq. ([4.13|) we now have the equation 

B 12 + B 13 + B 23 ]^(r 12 , r 13 , r 23 , R) = E^°\R) + £«(r 12 , r 13 , r 23 , R), (4.23) 

where now the function -E(r 12 , r 13 , r 23 , R) is given by equation ( |4. 1 0| ) but with the 
sum on B containing now one less operator. We seek solutions that are symmetric 
with respect to permuting the pairs 12, 13, 23. The leading inhomogeneous 
solution can be found as before, cf.( |4.1(j| ), 



rfUnz, rxs, r 23 , R) ~ S A (R) ^ + J + ^ . (4.24) 
The homogeneous solutions will be of the form ( 4.20 ), but with additional sum- 



mations over I" and m" for the third factor gi» m " as a function of the third small 
vector distance. This is all the information needed for the calculations that follow. 



C. Calculation of K(R) 

According to our strategy we need now to evaluate the limits shown in 
Eq. (|4.1j ). We can firstly compute the derivatives operating on the products of 
T 2 and discover that the limits are not singular. In fact that contributions decays 
rapidly in R like R~ 2C - h . In computing the derivatives of the cumulant of JF 4 we 
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need to consider all the contributions to the solution of JF 4 for two pairs of coalesc- 
ing points that were examined in subsection [4 A 2j , and seek the leading one. The 



solution (f4.15|) does not survive the derivative in Eq.(fO). The I — I' — term 



in (14.201) is constant and also does not survive. The next term that survives the 



derivatives is the / = /' = 2 term, and it is less singular than the inhomogeneous 
solutions ( |4.18| ) since /3 2 > (2- Accordingly the leading contribution is (|4.18|) , and 
upon performing the derivatives we get 

^)~%P lim t^f- ( 4 - 25 ) 

where we have used the fact that (h = 2 — ( 2 - The singular limit has to be 
understood in light of the full equation for T%, Eq.( |2.59|) , in which the K-diflusive 
terms are explicit. The role of these terms is precisely to truncate the divergence 
that is implied by ( f4.25| ) . As a consequence the divergence is only applicable in the 



inertial range with r 12 ,r 34 > rj, whereas in the dissipative regime the divergence 
disappears. Thus for evaluating K(R) via inertial range values we must replace 
the limit ri 2 , r 34 — > by r i2 = r 34 = r\: 

*W ~ • < 4 - 26 > 

Comparing to Eq. (|1.16|) we see that we recover the result that A = A c = ( h . 
We can rewrite ( }4.26|) in a final form by introducing e as 

e = k(|VT| 2 ) = -Kkm ViV 2 F(ri 2 ) oc n/rj Ch , (4.27) 

where we used the fact that ^ r (r 12 ) ~ r 4 |. Using the last equation in the preceding 
one we find the final result: 

K(R) ~ e 2 S 4 (R)/S 2 (R) 2 . (4.28) 
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D. Calculation of the correlation functions L;;/ m 

In this subsection we turn to the calculation of correlation functions that 
expose the scaling properties of anomalous fields with other irreducible repre- 
sentations of the rotation group. We can do this by introducing the following 
correlation functions: 

Li,i>,m(r 12 ,r u ,R) = J d cos 9 12 d cos 34 d(0i2 - 034) 
x .F 4 (r 12 , r 34 , R)Y ltrn (9 12 , MY v< _ m {6 M , 34 ) • (4.29) 

For I = V = the leading contribution to this correlation function arises from 
ACTED , and 

£0,0,0(7-12, r 34 , R) ^ ffL « # C4 " 2Ca • (4-30) 



For 1,1' > 2 the leading contribution arises from ( |4.2(J| ). Using ( 4.21|) we write the 
final result 

L lAm (r 12 ,r 34 ,R) ~ 5 4 (i2) (^)* • (4.31) 

Finally, for I = and /' > 2 or vice versa the leading contribution stems from the 
solution ( |4.22| ). The expression for A),z',o( r i2, r 34 , R) can be obtained from ( |4.31| ) 
simply by replacing (3 = by £2- 

Eq.( [4.3lD gives rise to a set of anomalous local fields. Taking gradients 
Vi Q V2/3, one produces one scalar field (Vi • V 2 ), which corresponds to I = and 
has anomalous scaling £ A = 2 — £ 2 , and the traceless tensor (Vi Q V 2 /3 — |ViV 2 ) 
corresponding to I = 2. Taking four gradients one can produce anomalous fields 
originating from I — 4, and so on. 



E. Calculation of the cumulant J%(R) 



The quantity J 4 (-R) has reducible contributions which were computed above, 
and a cumulant part which is obtained from JF£. The calculation of J4 (R) follows 
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very much the lines of the calculation of K(R), except that one needs to find again 
which of the solutions found in sec [4 A 2| contributes mostly to equations Q4.5|) and 



It turns out that the leading contribution to (|4.q ) stems from the solution 
( |4.15| ), whereas the leading contribution to (|4.6|) arises from the solution ( [4.241 ). 



As before the derivatives with respect to r x and r 2 result in a singular limit when 
7*12 — > 0, and we have to take T\ 2 = V after computing the derivatives. On the 
other hand the other limits are regular and they do not require special care. The 
evaluation of J^i(R) and J^ 2 {R) turn out to be the same up to constants. The 
result is 

Jt{R) ~ KSi(R)/R^rf h . (4.32) 

Using Eq. (|4.27p this can be written finally as 

Jt{R) = C 4 eS 4 {R)/S 2 {R) , (4.33) 

where C4 is a dimensionless constant that will be determined below. Together 
with the reducible contributions that were computed above we can write 

J 4 (R) = eS 2 (R) + |[V5 2 (i?)] 2 + C A eS 4 (R)/S 2 (R) . (4.34) 

It is obvious that the second term is small compared with the first and it can be 
neglected. 

5. THE CALCULATION OF J 2 n(R) 



The correlations ( 1.10|) which make J 2n {R) have Gaussian decompositions and 
a cumulant part. The leading contribution to the Gaussian decomposition is 



J2n(R) = ^S 2n - 2 (R) 



(5.1) 



33 



which corresponds to the first term in ( f4.34|) . To evaluate the cumulant part one 



needs to compute correlation functions of the type 

J P>q (R) = ((iviecroi^croe^n))) - (5.2) 

with i"i = — R/2 and r 2 = R/2. The calculation of J PA {R) follows from the 
equation for T 2n upon coalescing a group of p points (denoted below as the a 
group) into the position —R/2 and a group of q points (denoted as the j3 group) 
into the position R/2. We start with Eq. (|2.57|) which in the inertial interval takes 
on the form 



2n p 2n 



I £ B aa , + £ *V + £ E B a p]^ n (r l ,...,r 2n ) = RHS . (5.3) 

a>a'=l f3>@'=p+l a=l /3=p+l 

As before the effective equation is obtained by grouping together a quantity 
E{{r aa '},{rpp>},R) = RHS - £ £ B a(3 J r ^ n (v 1 , . . . , r 2 „). To find the form of 
the solution we note that when we had one pair of coalescing points this led 



to the solution (|3.14j). Two pairs of coalescing points led to (|4 . 1 6|) , whereas 



three pairs of coalescing points resulted in ( f4.24p . In the present case we have 



[p(p — 1) + q{q — 1)]/2 coalescing pairs and the solution which belongs to the same 
family is written as 



2n 



F2n^~S 2n {R)[ £ rt,+ J2 r%]/R^ . (5.4) 

a>a'=l P>/3'=p+l 

Next we need to compute the derivative with respect to r x and r 2 and take the 
limit of all pairs of coalescing distances going to zero. The only divergence will 
be associated with h whereas all the other limits are trivial. According to 
our strategy we have to cut the divergence at r 12 = rj and thus we calculate 
^2n,inh ~ S2n{R) / rf h ■ Using again Eq.( |1.17|) we have 



J c 2n (R) = C 2n eS 2n (R)/S 2 (R) , (5.5) 

where C 2n is an unknown dimensionless coefficient that will be determined soon. 
This equation is the generalization of ( ^4.33| ) for any n. Combining now the cu- 
mulant with the leading reducible contribution we get the final result 
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J2n(R) = £ S 2n _ 2 (R) + C 2n S 2n (R) / S 2 (R) 



(5.6) 



6. SCALING EXPONENTS OF THE STRUCTURE FUNCTIONS 

In this section we collect all the results obtained above with the aim of reaching 
conclusions about the scaling exponents of the structure functions. The first result 
that we need to pay attention to is ( ^4.26| ) or ( |4.28| ) for K(R). This result shows 
that K(R) ~ i?^ 4-2 ^ 2 . Since K(R) cannot be an increasing function of R we 
conclude immediately that 

C 4 - 2 < 2( 2 • (6.1) 



Next we need to consider (|4.34|) for J^(R), rewriting it after neglecting the second 
term as 

J 4 (R)^eS 2 (R)[l + C 4 ^\. (6.2) 

The second term in the parenthesis is dimensionless, and may be written as 
(£/R) 2 ^-<*. 

Now there are two possibilities: 

(i) C4 = 2C2 and the scaling is normal. The first and third terms have the 
same scaling and they are of the same order. 

(ii) £ 4 < 2( 2 , and the scaling is anomalous. If so, the second term in ( |6.2|) 
must be larger than the first. For that to happen the renormalization scale must 
be the outer scale L. 

The implications of the first possibility are somewhat strange. For example 
if there is normal scaling the correlation function K(R) does not decay in R. 
This means that the dissipation field is not mixing. On the contrary, if there 
is anomalous scaling, then the correlation K(R) decays as is expected from a 
random field. We will explore now the second possibility and show that if there 
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is anomalous scaling then the scaling exponents can be computed in agreement 
with Kraichnan's arguments. 

If we accept anomalous scaling, then the leading term in ( p.6| ) is the second 
one. Generally speaking J2 n (R) may be written (following Kraichnan et al. [0) 



as 



J 2n {R) = Ann J dA0P(A0)[A0] 2n ~ 2 (|V0| 2 |A0 



(6.3) 



where (|V0| 2 |A0) is the average of | V0| 2 conditional on a given value of A0. 
P(A0) is the probability to observe a given value of A0. Eq. ( |6.3|) is exact. The 
question now is what is the dependence of (| V0| 2 |A0) on A0. In order to recover 



our result (5J5) this conditional average must satisfy 



V0| 2 |A0) = C(AQ) 2 /S 2 (R). 



(6.4) 

This means that the coefficient C 2n in ( |5.5| ) is n- independent. We can determine 
this coefficient from the particular case n = 1. Since J 2 = e the conclusion is that 
C n = 1. Equipped with this we consider again the balance equation ( |1.8|) 



Rl ~ d M Rd ~ lh{R) M S2n{R) = j2n{R) 



(6.5) 



Using now the definition of the scaling exponents S2 n {R) ~ R^ 2 " we retrieve from 



the result ( |1.12| ) which can be also written as 

(2n=\ (2~d+ V /(C2 + rf) 2 + 4C 2 (^-l) 

This is Kraichnan's anomalous scaling. 



(6.6) 



7. SUMMARY AND CONCLUSIONS 

The first conclusion of this paper is that renormalized perturbation theory 
for hydrodynamic fields has the potential to describe nonperturbative effects. 
Exact resummations of the diagrammatic series result in exact equations for the 
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statistical quantities that contain not only perturbative but also nonperturbative 
effects. 

Secondly, the passive scalar problem with rapidly decorrelating velocity field 
displays a particularly simple analytic structure in which all the scaling properties 
emanate from one differential operator B(R). As a result this theory becomes 
"critical" in the sense that the anomalous exponent A that was introduced in 
paper I is exactly critical. The reason for this is that both A and A c come from 
the same operator, and thus they must be the same. In such a situation the 
subcritical scenario that was suggested in || is untenable in this case. 

Thirdly, the simplicity of the theory allowed us to compute the whole spectrum 
of anomalous exponents that are associated with the ultraviolet divergence. It 
was shown that these exponents are related to the spectrum f3i ( |3.11| ) which 
appears in the law of isotropization of the 2-point correlation function. The same 
exponents determine the R dependence of the correlation functions Lw >m that 
were introduced in ( |4.29| ). 

Finally we showed that as far as the structure functions S n (R) are concerned, 
there are only two possibilities. Either the exponents satisfy <^2n = n^2 and the 
scaling is normal or the scaling is anomalous with the outer renormalization scale 
and with the law (|6.6|) . We noted that normal scaling also implies that the 
dissipation field is not mixing. 

Whether or not the "critical" situation of this model with a marginal scaling 
exponent A = A c is structurally stable is an extremely important question that 
must await further research. 
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APPENDIX A: DERIVATION OF THE EQUATION FOR THE SIMULTANEOUS 

TWO-POINT CORRELATOR 



Here we obtain the equation of motion for the simultaneous two-point corre- 
lator, .F(ri, r 2 , t — 0). We make use of Eq. ( |2.18| ) to verify that the operator 



G 12 (r |ri,r 2 ,t) = d t + V 1 (r 1 - r ) + V 1 (r 2 - r ) 



is the inverse of the product of Green's functions, that is, 



(Al) 



Gu( T o\ri> r 2> t)g(r \r u r[, *)</(r |r 2 , r 2 , t) 
= S{T 1 -T f 1 )6(r i -i f 2 )6(t). 



(A2) 



We have used here the fact that the Green's function when multiplied by S(t) 
may be evaluated at t — 0. Applying this operator and performing the spatial 
integrations, one obtains directly 



Pi(ri - r ) + X>(r 2 - r ) T{t u r 2 ) = $(r |ri, r 2 ) + $o(ri - r 2 ). 



(A3) 



Now using the definition of $(ro|ri, r 2 ) from Eq.( |2.12 ) and dropping the ro de- 
pendence, one may rewrite this as 



8 3 



where the operator H is given by 

/ d 
H(t 1 ,t 2 ) = \hi j (r 1 )-^—+h ij (r 2 ) 







^"(ri,r 2 ) = $ (ri,r 2 ) 



d d 



(A4) 



(A5) 



dr 2i ) \drij dr 2j J 

In the case that all quantities are only functions of R = ri — r 2 , ft vanishes. 
The remaining terms are equivalent to the definition of the operator P 2 (R) in 
Eq. (|2.32j ) and one recovers Eq. (|2.29| ). 
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APPENDIX B: EQUATION FOR THE CUMULANT OF THE 2iV-POINT 

CORRELATOR 



In this appendix we derive equations for the many-time correlators using the 
diagrammatic approach. The infinite diagrammatic series for the 4-point corre- 
lator was presented in M . Let us consider the series for the correlator. A typical 
diagram for a 2n-point correlator consists of n 'tramways': strings of Green's func- 
tions, connected in pairs at one end by two-point correlators. These tramways 
may be interconnected by velocity correlators. Therefore to build up an arbitrar- 
ily complex diagram, one adds successively connections between chosen pairs of 
tramways. Let us consider a given pair of end-points, x a and xp. Group all dia- 
grams together in which the Green's function beginning at x a and that beginning 
at Xp are linked by a velocity correlator. The series of diagrams to the right of 
this first velocity correlator is again the series of diagrams for the full correlator. 

The difficulty with writing a resummed equation for the full correlator is that 
the lowest order terms are Gaussian, and disconnected. In building upon the 
Gaussian terms by the addition of 'rungs' as just described, one will generate, 
amongst other terms, a series in which the disconnected parts remain discon- 
nected, and may be resummed again into the Gaussian decomposition. This 
means that these terms actually appear twice. In order to avoid this one may 
write an an equation for the cumulants. For the cumulant the lowest order terms 
are those in which all tramways have one velocity correlator connection to an- 
other. 

To simplify the appearance of this equation we introduce the opera- 
tor C a p which represents the addition of a rung, and here operates on 
•7"2n(0|xi, x a , Xp, X2n)- The definition is 

C a p * J 7 2n(Q\Xl,:.,X a ,Xp,...,X 2 n) 
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dr[dr' 2 / dtQl{Q\x a ,xp,r' x ,t,r' 2 ,t) 



(Bl) 



X Hij(r[, r '2)~Q^F~Q^F~F 2n(0| a; l? •••) r 'u r 2' X 2n)- 



Since JF 2n (0|xi, x 2 , .., x 2 „) is symmetric with respect to all exchanges of coordi- 
nates this definition is sufficient for any pair of indices 1 < a, (3 < 2n. 

The lowest order terms may be expressed in terms of the operator C and 
two-point correlators as 

2n 2n-l 



^2n,oi^\ x U x 2, -, x 2n) = II C M+1* II F{ x k, x k+l) 

j=2,4 fc=l,3.. 

Then the resummed equation has the form 



(B2) 



perm 

+ ^2C a p*J 7 2n {0\x 1 ,X2,..,x 2n )- 

a>f3 



(B3) 



Now operate on both sides of the equation with the product of the inverse Green's 
functions Q2.26 ). Using the fact that 



^- + Airaij \Jj- + Airpfj C al3 = 5(t a - tp)B a p, 



(B4) 



one finds 



2n ( d \ 

V2n-2 



£ £ 

<«,/?> (7i) 

perm(*yi ^a,p) 



n 6 ( l h -hi+i, 



i=l 



_d_ 

dt a 



' d_ 



+ t>i(r (3 ) B 7l)12 J 7 (x a ,x 



2n-4 



^7j+i.7j+2^ r ( r 7j' r 7j+i) ( ^7j+l7j+2*^ r ( a '72n-2 3 x p) 
,5=2,4.. J 

) x 2 j •■ 1 x 2n ) ■ 

a>/3 7^a,/3 



(B5) 



Therefore we have derived a closed equation for the time development of the 
many-time cumulants of the 2n-point correlator in terms of the operator B and 
the two-point moments only. We will not solve this equation in this paper. 
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